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^ : Abstract 

In a recent experiment the superconducting gap of a single unit cell thick FeSe film on SrTiOa 
substrate is observed by scanning tunneling spectroscopy and angle-resolved photoemission spec- 
troscopy. The value of the superconducting gap is much larger than that of the bulk FeSe under 
ambient pressure. In this paper we study the effects of screening due to the ferroelectric phonons 
on Cooper pairing. We conclude it can significantly enhance the energy scale of Cooper pairing 
and even change the pairing symmetry. Our results also raise some concerns on whether phonons 
can be completely ignored for bulk iron-based superconductors. 
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I. INTRODUCTION 



In a recent experiment [1] a single unit cell thick FeSe film is grown on the Ti02 terminated 

(001) surface of SrTi03(ST0) by molecular beam epitaxy. Two gaps (~ 10 meV and 20 meV) 
are observed by scanning tunneling microscopy (STM) at low temperatures. At present there 
is no transport data and Tc is not determined. Subsequently an angle-resolved photoemission 
spectroscopy (ARPES) result appears. [2] According to this report, there is only electron 
Fermi surface, suggesting FeSe is electron doped. The observed electron pockets are nearly 
circular, and an approximately constant (~ 15 meV) superconducting gap is detected on 
them. Significantly, by studying the temperature dependence of the energy gap, an estimate 
of Tc, namely, 55 ± 5 K, is obtained. Currently there is no explanation of the two gap versus 
one gap discrepancy between the two measurements. 

There are many unanswered questions. For examples, (1) What caused the FeSe doping? 

(2) Can the discrepancy between Ref.[l] and Ref.[2] be due to surface doping caused by 
the sample treatment prior to the ARPES measurement? (3) Are there buried, hence not 
yet detected, interface metallic bands? [2] (4) Are there interface ferroelectric ordering? (5) 
How strong is the coupling between the electrons in FeSe and the ferroelectric (FE) phonons 
in STO? Due to the uncertainty of the surface doping we shall study both undoped and 
electron doped FeSe. 

STO is a "quantum paraelectric" insulator. The huge dielectric constant at low temper- 
atures is due to the ionic movements. The associated phonon, the FE phonon, involve the 
relative displacement of the Ti and O atoms. An example of such displacement is shown in 
the right panel of Fig. 1. This phonon mode is soft at the zone center. At low temperatures 
bulk STO also undergoes an antiferrodistortive (AFD) transition. This distortion involves 
alternating clockwise and counterclockwise rotation of oxygens about titanium (see the left 
panel of Fig. 1) In a recent local density functional study [3] it is found, aside from zone 
folding, even static AFD distortion has little effect on the bandstructure. Therefore we shall 
ignore the AFD phonon in the following discussions. 

In the rest of the paper we perform a two-stage renormalization group (RG) study. [4, 5] 
The first stage is a functional renormalization group (FRG)[6, 7] calculation. It determines 
the most important electronic scattering processes at low energies. As found in previous 
studies [7-9] for energies lower than the magnetic fluctuation energy scale, Ag, the strongest 
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FIG. 1: Left panel: an example of AFD distortion in the Ti02 plane. Right panel: a caricature of 
the frozen FE phonon with the atomic displacement in the (001) direction. 

electron-electron scattering is in the Cooper channel (for superconducting samples). This is 
used as the input for the the second stage analysis. Here we study the effect of FE phonons 
at energies lower than the maximum phonon frequency Ap/j < Ag. This is done by solving a 
generalized Eliashberg equation. The FE phonons have two important effects on the FeSe 
electrons. The first is to mix states separated by momentum (vr, vr) in the unfolded Brillouin 
zone. This is because a frozen FE phonon breaks the z o —z reflection symmetry, the 
symmetry that enables the Brillouin zone unfolding. In addition, the soft zone center FE 
phonons can screen the intra-pocket electron-electron repulsion. 



II. THE MAIN RESULTS 



To present the underlying physics, let's simplify the electronic structure by retaining 
only two Fermi pockets: the hole and electron pockets for the undoped FeSe/STO, and the 
two electron pockets for the doped FeSe/STO. The two main effects of FE phonon discussed 
earlier is described by an intra-pocket and an inter-pocket coupling constant, Aintra and Aintcr- 
The superconducting transition temperature determined from the two-stage RG approach 
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is given as follows. 



Odd sign pairing : Tc = Aph ■ exp 



Even sign pairing : Tc = Aph ■ exp 
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where 



A± Ajntra i Ajntor 
V± = Mntra i Venter 

K^ = \4/[l + \41n(Ae/Aph)]. 



(2) 



Here "Odd" /"Even" sign pairing means the gap function has opposite/same sign on the two 
Fermi pockets, ^intra and Venter are the average of the intra and inter-pocket Cooper scattering 
strength obtained from the first stage calculation. Xmtra and Xmter are the strength of the 
phonon mediated Cooper scattering entering the second stage calculation, (u) is a weighted 
average[10] of the phonon frequency satisfying < (u;)/Aph < 1. In obtaining the above 
result we have assumed that 1 + V±ln(Ae/Aph) > 0, i.e., the pure electronic driven Tc os 
lower than Aph/fc^. 

According to Eq.(l) increases Ajntra — Aintcr, and Venter — ^intra raises the for the odd 
sign pairing. On the other hand, to raise the Tc of the even sign pairing we need to increase 

Aintra + Ainter but deCreaSC Vintra + V'intcr, l-G., 



(Aintra ~ Ainter) and (Venter ~ ^iitra) 

(Aintra + Ainter) and (Venter + ^intra) \( 



rpodd 

rpcvcn 



(3) 



The physics behind Eq.(3) is rather simple. Odd/Even sign pairing requires the inter-pocket 
Cooper scattering to be repulsive/attractive. Since phonon mediated scattering is necessarily 
attractive, it follows that strong inter-pocket electron-phonon interaction enhances even 
sign pairing while suppress the odd sign pairing. Because attractive intra-pocket scattering 
strengthen both even and odd sign pairing, strong intra-pocket electron phonon interaction 
is beneficial to both. 
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FIG. 2: (a,b) the phase diagram. Light blue/gray denote even/odd sign pairing, respectively. (c,d) 
the Tc enhancement factor Tc/Tcq where Tco is the superconducting transition temperature in the 
absence of the electron-phonon interaction. The parameters we used to construct the figures are 
Ae/Aph = 2, {u})/Aph = 0.5, and Tenter = 0.2, T^ntra = 0.05 for (a) and (c), and Tenter = 0.2, 
Mntra = —0.05, for (b) and (d). The small triangular region near the lower left hand corner of 
panel (a) is non-superconducting. 

A. The undoped FeSe/STO 

For undoped FeSe/STO the first stage RG generates Venter > hence favoring odd sign 
(S'+_ ) pairing. [7, 8, 11] This is caused by the antiferromagnetic fiuctuation. As discussed 
above, Aintra enhances the pairing while Ainter weakens it. Setting Venter = 0.2 and 
Vintra = ±0.05, the "phase diagrams" as a function of Aintra and Ainter is shown in Fig. 2(a,b). 
The associated Tc enhancement factor is shown in Fig. 2(c,d). Note the magnitude of the 
Tc enhancement differs by approximately an order of magnitude in Fig. 2(c,d) by merely 
reversing the sign of Mntra- Clearly this quantity is not something we can confidently predict. 
What is robust is the fact that when Ainter >> Aintra electron-phonon interaction stabilizes 
even sign pairing. Conversely for Aintra >> Ainter odd sign pairing is favored. Since for 
undoped FeSe/STO the FE phonons mainly cause intra-pocket electron scattering ( Aintra >> 
Ainter), wc cxpcct the electrou-phonou interaction to strengthen the odd sign, in this case 
S^^, pairing. 
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FIG. 3: Left: A tight-binding fit to the band dispersion observed in Ref.[2]. Right panel: the 
electron-Uke Fermi surfaces. 

B. The electron doped FeSe/STO 

Our tight-binding fit of the ARPES bandstructure and the associated Fermi surfaces are 
shown in Fig.3. Because the hole bands have completely sunk below the Fermi energy, the 
antiferromagnetic fluctuation only occurs at energies greater than the separation between 
the top of the hole bands and the Fermi energy. According Ref. [2] this separation is ap- 
proximately 80 meV. If sufficiently strong this high energy magnetic fluctuation can trigger 
5'++ pairing on the electron pockets. This pairing form factor can be thought of as 5*+- re- 
stricted to the exposed Fermi surfaces. When this form factor is the leading pairing channel, 
we expect Mntra and Venter to be both negative. Under this condition adding the electron- 
phonon interaction the phase diagram is shown in Fig. 4(a)- there is only 5*++ phase. The 
Tc enhancement is shown in Fig. 4(b). 

For larger separation between the hole bands and Ep, the nodeless d-wave (where the 
gap function has opposite sign on the two electron pockets) is the leading pairing channel. 
As pointed out in Ref. [12] in the presence of hybridization between the electron pockets (due 
to the absence of the z ^ —z symmetry at the interface) the d-wave pairing can become 
nodal. For strong hybridization the reconstructed electron pockets could have opposite sign 
pairing due to the repulsive electron-electron interaction. [12, 13] We expect in both cases 
Fig. 2 should apply, i.e., strong enough inter-pocket scattering (the first main effect of the 
FE phonon in earlier discussions), can destabilize the odd sign pairing and turn it into even 
sign. In addition to the effect of phonon the inevitable disorder scattering also tends to 
destabilize the odd sign pairing in favor of even sign. Thus we strongly believe the electron 
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FIG. 4: (a) the phase diagram. Light blue denotes the even sign pairing, (b) the Tc enhancement 
factor TcjTffy where T^a is the superconducting transition temperature in the absence of the electron- 
phonon interaction. The parameters we used to construct the figures are Ag/Ap/i = 2, ((x!)/Aph = 

0.5, and Mntra = -0.2, Mnter = "0.05. 

do-ped FeSe/STO have even sign (or S^^) pairing. 

Now we present the more technical part of the paper. 




III. THE FIRST STAGE - ERG CALCULATION 



The electronic Hamiltonian consists of the bandstructure part, -ffband, and the Hubbard- 
Hundes interaction Hi. The details of Hj can be found in Appendix A. For the undoped 
FeSe/STO because of the well nested Fermi surfaces [3] particle- hole scattering is expected to 
grow continuously at low energies. As the result the assumption made earlier, that we only 
need to retain the Cooper channel electron-electron scattering for E < Aph, is not a priori 
justified. Consequently in the next subsection we shall perform a full FRG treatment of 
both electrons and FE phonons to check the results obtained using the Eliashberg approach. 



A. The undoped FeSe/STO 

The Fermi surfaces (in the unfolded Brillouin zone) of the undoped system are shown in 
Fig. 5(a). The phonon Hamiltonian is given by iJph = X^p '^(p)'^p'^p- Here p is the three- 
dimensional momentum. We assume the following dispersion for the FE phonon a;(p) = 
a/wq + c^p^ where Uq ~ 2meV and c ~ 70meV ■ A are estimated from an early neutron 
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FIG. 5: (a) FeSe Fermi surface(s). (b) The superconducting critical scale Ac versus Amtra- (c) and 
(d) presents the flow of Ssc and Ssdw, respectively. (Higher scales are not shown). The value of 
Aintra increases by equal intervals along the arrows. 

measurement. [14] The electron-phonon Hamiltonian is given by 

i,a Pz ka,k'f) 

Hc-ph describes the coupling between electrons in the FeSe layer and the nearest ion dis- 
placement Ui in the Ti02 layer of STO. In the above = pzZ is the out-of-plane momentum 
of the phonon and is the third Pauli matrix in the Nambu space. The summation over 
follows from the fact that the elelctron-phonon coupling occurs at the interface. Here a and 
b label the electron bands, and 5'ka,k'b;p2 where (ka|k'6) is the 

overlap between band Bloch states. The phonon mediated intra and inter-pocket Cooper 
scattering strength are given by Aintra = Y1=i and Aintcr = where 

where ((■)) denotes the joint average over k and k' which lie on Fermi pockets a and 
b respectively. Using uq and c above and the Fermi pockets in Fig. 5(a) we estimate 

AinterZ-^intra ~ 1/12. 

We generalize the SM-FRG method[15, 16] to treat the effects of both electron-electron 
and electron phonon interactions. In principle, we envision a boson-fermion FRG calculation 
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that involves the flow of electron self energy, the phonon self energy, the elelctron-phonon 
coupling vertex, the 4-point phonon vertex, as well as the electron-electron interaction ver- 
tex. However in the following shall view the electron dispersion, the phonon dispersion, 
and the elelctron-phonon coupling as fully renormalized quantities. The flrst two can be 
determined from experiments and the last quantity is viewed as an adjustable parameter in 
our theory. What's left is the RG flow of the electron-electron interaction including con- 
tributions from the pure electron-electron interaction and the phonon mediated interaction. 
In the following the renormalized interactions in the superconducting (SC) and spin density 
wave (SDW) channels are denoted as Ssc,sdw- (The charge density wave channel turns out 
to be unimportant and will not be discussed.) The deflnitions as well as technical details of 
the SM-FRG method can be found in Appendix A. 

Fig. 5(b) summarizes the superconducting critical scale Ac (flUed squares) versus the 
electron-phonon interaction parameter Aintra (Amtcr ~ Aintcr/12). This is extracted from 
Fig. 5(c) which shows the flow of Ssc- Fig. 5(b) shows a linearly rising Ac as a function 
of Aintra- When we compare the maximum Ac with that in the absence of electron-phonon 
interaction a maximum Tc enhancement ~ 6.5 is obtained. Given the fact that Aintra ~ 
12Aintcr this result is consistent with Fig. 2. We have checked that for all values of Aintra 
in Fig. 5(c) the pairing symmetry remains to be S+-. Fig. 5(d) shows that the Ssdw 
almost saturates at low energy scales, and the electron-phonon coupling reduces 1 / Ssdw only 
slightly. Since this type of SDW is related to the pairing interaction, the above observation 
justifles our previous assumptions for the input to the Eliashberg calculation. 

B. The electron doped FeSe/STO 

In this case the assumption of only retaining the Cooper channel electron-electron scat- 
tering is a good approximation. Hence we shall focus on the pure electronic calculation to 
flnd out the leading pairing channels and use it as input for the second stage Eliashberg 
calculation. The electronic structure we use to describe the electron doped FeSe is shown 
in Fig. 3. It it worthy to note that this bandstructure differs substantially from those for 
other pnictides in that the hole band top near F and the electron band bottom near X and 
Y are separated by a small gap. Interesting in a recent ARPES work on AxFeySe2[n] this 
feature is noted and emphasized. 
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FIG. 6: Left panel: the RG flow of the Cooper scattering strength Ssc associated with the two 
leading pairing channels. Black circles: s-wave pairing. Red squares: dxy pairing. Mid panel: the 
s-wave form factor. The polar angle 6 is defined with respect to the +x direction for both X{7r,0) 
and Y{0, vr) electron pockets. Right panel: the d^y form factor. 

For {U,U',Jh) = (2,1.7,0.15) eV, the FRG results for the Cooper pairing is shown in 
Fig. 6. The left panel shows the leading pairing instability being in the s-wave channel 
(black), while the subleading pairing channel has dxy symmetry (red). The form factors 
for the s and d^y pairing are shown in mid and right panels respectively. For the s-wave 
pairing, the form factor has the same sign on the two electron pockets X and Y. However 
the amplitude of the form factor is quite anisotropic. In contrast the dxy form factor has 
nodes on the two electron pockets. We have also checked for considerably weaker interaction 
parameters, the leading pairing channel is dx^-y^ where the gap function has opposite sign 
on the X and Y electron pockets. As discussed earlier, the interaction with FE phonon tends 
to stabilize the s wave pairing and suppress the d-wave pairing. For intermediate electron- 
electron interaction parameters we obtain s-wave pairing. The addition of FE phonons will 
further enhance the pairing energy scale. 



IV. THE SECOND STAGE - THE ELIAHBERG EQUATIONS 

The formalism and notations in this section closely follow those in Scalapino et a/. [18] 
and McMillan [10]. The details can be found in Appendix B. In order to make analytic 
calculation feasible, we model the low energy electronic degrees of freedom by two Fermi 
surfaces. In the case of undoped FeSe these correspond to hole and electron pockets and 
for the electron doped case both Fermi surfaces correspond to electron pockets. These two 
Fermi surfaces are labeled by a = 1, 2 in (Bl) respectively. Moreover we shall assume these 
Fermi surfaces have constant density of states(DOS) A'^i^2- In addition, in the rest of this 
section the interaction Vabcd{Pi, Ps^ Pd will be assumed to depend on the Fermi surface 
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index only, i.e., Vahcd- These simplifications are made in view of the presumably strong 
disorder-induced quasiparticle scattering at the FeSe/STO interface which will average out 
the fine structures in Vabcd{Pi^P2iP'i^PA) ^'^^ ^^e DOS. 

Assume the following form for the self energy of the Nambu spinor ^, 

T.a{uj) = [1 - Za{uj)]uTQ + A„(a;)ri, (6) 

we derive the self-consistent equations for Z and A following standard procedures. [18] The 
results are given by Eqs. (B7) and (B8) in Appendix B. Assuming the McMillan ansatz 
(Eqs (B9)) these equation can be solved to yield ^^(0) = 1 + X];, a/ Nb/NaXab for the normal 
state Za{0), and the following eigenvalue problem for Tc(B8) 



6=1 



In 



A 



ph 



ah 



A 



ph 



2 



In- 



A 



ph 



^Af,(0) + In 



A. 



Aph 



iVhAJoo) 



A^bAb(cx)) 



/iV^A,(cx)) = v^Z,(0)A„(0), 



(7) 



'NaK{00). 



In the above equations \ab = 2 and where F{u) is the density of states 

associated with the phonon mode and a^^bi^) is the effective electron-phonon coupling 
constants. In addition Vab = V ^aNbVabba, and {u)ab is a weighted average [weighted by 
Q!^j(i^)F(z/)/i^] of the phonon frequency. See Appendix B for details. 



For the sake of simplicity we shall set A'^i = N2, Xaa A 



intra 



(Aii + A22)/2, Xab — ^ Ainter 



(A12 + A2l)/2, Vaa ^^intra = (^^11 + ^^22)/2, Vab "> ^^intcr = (^^12 + ^^2l)/2, and {uj) ab = (c^) in 

the following. With the above simplifications we can determine the Tc for 5*+- (Ai = — A2) 
and 5*++ (Ai = A2) pairing as given by Eq. (1). We have checked changing N1/N2 has little 
effect on while have a strong effect on A1/A2. 

If phonons are indeed responsible for the observed enhancement, there should be other 
signatures of the electron-phonon coupling. For example due to the large anharmonicity of 
the FE phonons we expect using a pump laser to excite them will have a clear effect on 
Tc- The more traditional phonon signatures such as kink in the quasiparticle dispersion and 
shoulder in the tunneling experiment are discussed in Appendix C. 



V. DISCUSSION 



We have studied the screening effects the ferroelectric phonons of SrTiOa on the inter- 
action between the electrons in FeSe. We conclude such coupling can enhance the pairing 
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strength of FeSe. Moreover we find when the inter-pocket electron-phonon scattering is 
strong, the opposite sign pairing will give way to the equal sign pairing. 

An immediate question one might ask is why don't the FE phonons have similar effect 
on the Tc of doped STO. We believe the answer is polaron formation - for a range of strong 
electron-phonon coupling, the formation of polarons instead of Cooper pairs are favored. 
In a recent optical experiment on n-type doped STO, [19] a very sharp Drude peak with a 
substantial mass enhancement (consistent with that of "large polarons") was observed. 

The current study raises the concern about whether the role of phonon can be com- 
pletely ignored in bulk iron-based superconductors. [20] With appropriate interpretation of 
Aintra and Ainter, our rcsults can be used to address the phonon effects in bulk iron-based 
superconductors as well. 

Material wise there are other nearly ferroelectric perovskite materials, for example 
KTa03.[21] If FeSe films can be epitaxially grown on these materials similar enhance- 
ment should occur. Finally the results of Ref.[l] and the present paper suggest the ... 
FeSe/(STO)„/FeSe/(STO)„... superlattice is a promising artificial material with high Tc. 
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Appendix A: The SM-FRG method 

In this section we describe the SM-FRG method we used. The microscopic hamiltonian, 
which is valid for all cutoff < bandwidth, we use is H = i^band + Hj + _ffph + -ffc-ph- Here 
-f^band is the two dimeusioual bandstructure. For the undoped case i/band is a courtesy of 
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Z-Y Lu.[22] For the electron doped case it obtained by a fit to the ARPES result [2]. Hj 
describes the local electron-electron interaction given by 

i,a i,a>l3 i,a>l3,a,a' 

+ Jh ^ (V^l,a,t^l,a,iV'i,/3,iV^i,/3,t + h-C), (Al) 

where V'i,o,o- annihilates a spin a electron at site i in orbital a{a = Sz'^—r'^, xz, yz, x'^—y'^, xy), 
ni,a,a = '^iaa'^i,a,a and Tij^a = Xlo-^i.a.o"- Calculation of the main text we used intra- 

orbital repulsion U = 2eV, Hund's rule coupling Jh = 0.31el^, and inter-orbit repulsion 
U' = U - 2Jh. 

Fig. 7 (a) shows a generic 4-point vertex function ri234 which appears in the interaction 
V'iV'2('~ri234)^/'3^/'4. Here 1,2,3,4 represent momentum (or real space position) and orbital 
label. The spin a and r are conserved along fermion propagators and will be suppressed 
henceforth. Figs. 7(b)- (d) are rearrangements of (a) into the pairing (P), the crossing (C) 
and the direct (D) channels in such a way that a collective momentum q can be identified. 
The dependence on all other momenta and orbital labels is written as 

Klt-K-r>,P+n ^ J]/:,(k,a,/3)P^.(q)/.(p,5,7), 

mn 

rk+i,k,p+q ^ E /-(k, a, 7)C_(q)/„(p, S, (3), 

mn 

rg^p,p+q,k ^ E/-(k,a,5)D_(q)/„(p,7,/3). (A2) 

mn 

Here /m=(i,o)(k, a, /3) = hi(k)Aio{a, (3) is a composite form factor, where hi(k) is chosen 
from a set of orthonormal lattice harmonics, A^o is a matrix in the orbital basis. 

The decomposition Eq. (A2) for each channel would be exact if the form factor set is 
complete. In practice, however, a set of a few form factors is often sufficient to capture the 
symmetry of the order parameters associated with leading instabilities. [15, 16] In our case, 
the lattice harmonics are chosen as /i(k) = 1, cos kx^cos ky, 2 cos kx cos ky and 2 sin kx sin ky. 
They are all even since only singlet pairing is relevant in our case. The A^-matrices are 
chosen so that the combination (pa-^la, P)4>fB {(pa is the real atomic orbital function) is 
irreducible and transforms according to Aig, Big or under the point group. [23] (One may 
also use any bilinear 0q0^ to determine a matrix A4, but it is less transparent symmetry 
wise.) Moreover, the Al-matrix is normalized as TrA^^A^ = 1. If the total number of 
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FIG. 7: A generic 4-point vertex (a) is rearranged into the pairing (b), crossing (c) and direct (d) 
channels. Here k, q, p are momenta, a and r denote spins which are conserved during fermion 
propagation, and m,n denote the form factor (see the text for details). 

composite form factors is A^, then P, C and D are all x matrix functions of q. Note 
that the P, C and D channels are not orthogonal. The overlap between different channels are 
important for the growth of pairing interaction out of, e.g., the magnetic interaction. [7, 9, 16] 
In the following we denote Xk = KX as the projection of X into the K-channel via Eq. (A2). 

In the case with electron-phonon interaction we need the phonon-mediated electron- 
electron scattering. According to the H^-ph defined in the text, this vertex is given by, 



^(q, I'n) OC 



OC 



V2 + a;2(q) 



(A3) 



where q = {q^, Qy) and z/„ are the momentum and (Matsubara) frequency transfer in the 
electron-electron scattering, and the last proportionality holds to leading order in ujq/c. The 
Aintra.intcr discussed in the main text are both proportional to A^. Notice that this vertex is 
naturally in the D-channel and bears trivial form factors with h = 1 and Ai = 6ai3 since 
according to Eq. (4) the electron-phonon interaction is local in real space and diagonal in 
orbital basis. 

The partial flows of P, C and D are given by the one-particle-irreducible diagrams shown 
in Fig. 8. Here the dashed line denotes 4-point fermion vertex, wavy line the (surface) 
phonon-mediated vertex, and the dash-wavy line means that both types of vertices can 
enter. We write the partial flow equations as, in matrix form and for a collective momentum 
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FIG. 8: One-loop diagrams contributing to the flow of the the 4-point vertex function in the pairing 
channel (a), crossing channel (b), and direct channel (c)-(e). Here m,m'n,n' denote form factors, 
while the momentum, orbital and spin indices are left implicit. The dashed line represent 4-point 
fermion vertex, dash-wavy line means that both fermion vertex and phonon propagator can be 
inserted but separately. The diagrams are one-particle-irreducible with respect to both fermions 
and phonons, and the phonon line shares the loop frequency since external fermion fields are set 
at zero frequency as usual. The Matsubara frequency is continuous but subject to hard infrared 
cutoff at running scale A. 



— = iC + Vc)x',HiC + Vc) 



dP 

dA 
dC 

dA 
dD 

— = {C + Vc- D)x',,D + Dx',,{C + Vc-D). (A4) 

Here A is the running Matsubara frequency cutoff, Vk is the phonon mediated interaction 
projected into the K channel at cutoff A, and ^ire matrix-kernels with elements 



1 f (Pp 

{Xpp)mn = -— y ^^^/r„(p,a,/3)Ga^(p + q,2A)G^5(-p, -iA)/*(p,7,5) +(A->-A), 
1 f cPt) 

ix'ph)mn =-^ J j^fUP,a,P)G^^{p + q,iA)Gsp{p,iA)f*{p,-f,S) + (A ^ -A)(A5) 

where G is the bare fermion propagator in the obital basis, and the summation over orbitals 
is left implicit. 
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Clearly, because of the A dependence, the effect of phonon mediated interaction, V, is 
important only if A reaches the phonon band, above which the main contribution to the 
flow of the fermion interaction vertex is from pure electron-electron interaction. However 
the electronic excitations do modify phonon self energy and elelctron-phonon vertex even 
when the cutoff scale is above the phonon bandwidth. In the present work such effects are 
accounted for by using the experimentally measured phonon dispersion and the elelctron- 
phonon coupling constant. The flows in Eq. (A4) collect contributions from independent 
one-particle- irreducible diagrams for the total change dT, which need to be subsequently 
projected to the three channels. Therefore the full flow equations can be formally written 
as, 



for K = P,C and D. We used the fact that KdK = dK by definition. 

The functions P, C and D are related to the effective interactions as, Vgc = —P — Vp 
in SC-, Vsdw = C + Vb in SDW-, and Vcdw = C + Vc-2D in CDW-channels. We monitor 
the most negative singular values Ssc,sdw,cdw of such interactions (for all q) versus the 
running cutoff A. The most negative one among Ssc,sdw,cdw tells us which channel is 
becoming unstable. The associated eigen function dictates the symmetry and wave vector 
of the order parameter. 

Appendix B: Multiple band Eliashberg equation 

In this Appendix we briefly outline the derivation of multiple band Eliashberg eqautions 
for the electron-phonon Hamiltonian (Bl), following the single band case of Scalapino et 



The effective elelctron-phonon Hamiltonian we consider, at the energy cutoff Ag, is given 




(A6) 



al. [18] 



by 



2 




a=l p q 




(Bl) 



a,b p,p' 



+ 2^^ Yl V-bcdiPl,P2,P3,P4)i%,cn'i^p,a)i'^l,d^3'^p,l>)- 



a,b,c,dPi,P2,P3 
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Here Ae is much smaller than the bandwidth but larger than maximal phonon frequency 
Aph, \l/pa is the Nambu spinor for electron band a with dispersion epa, aq is the destruction 



operator for phonon with dispersion uj{q), 



+ P4 = Pi + P2 - Ps' ^3 is the Pauh 



matrix, and Vabcd{PiiP2iP-iiPA) is the (effective) electron-electron interaction. Different 
levels of simplification will be applied to Eq. (Bl), which will be discussed in more details 
in the following. 

Assume the Green's function of Nambu spinor ^ to be given by 



The self-consistent equation of the self-energy w) is 



(B2) 



E„(p, iu) = -T^^ rsGbip', iuj')T-i ■ \gpp',ah? D{p - p', iu - iu') 

u',p' b (^B3) 

+ Vabba{P,P' ,P\P) 

where u' is fermion Matsubara frequency. D{q^v) is the Green's function of the phonon. 
Use the spectral representation 



1 



1 



iv' — V iu' + v 



D{q,iu') = / duB{q,u) 
Jo 

and sum over u' by the procedure of Ref. [18], the self-consistent equation becomes 



(B4) 



X \gpp',ab? I duB{p - p', u) 



N{u) + f{-uj") , N{u) + f{uj") 



u — 00 — u 



OJ — OJ" + V 



(B5) 



1 r+oo 1 

Mr^Gb{p\ uj")n]-Vabba{p - p') tanh(/3a;72) 
b ^ p' ^ 

Assume each electron band a has a circular Fermi surface with constant DOS Na-, define 



o?fJyy)F{y) as the average of \/Nal^b\9ppi ,ab^B{jp — p', v) over p on Fermi surface a and p' 
on Fermi surface 6, and ignore the momentum dependence of E close to Fermi surface, this 
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equation further simplifies to 



TT 



X 



du" / duali^{u)F{u)lm[T3Gb{p\uj")T3,] 
N{u) + f{-u") , iV(z/) + /(a;")- 



(B6) 



a; — (jj — z/ 

-oo 



U — tU + U J 



- - / du"lm[T^Gt{p',oj")n]-Vabt8inh{/3u"/2) \. 

Here N^u) = l/{e^'^ — 1) and /(w) = l/{e^'^ + 1) are the Bose and Fermi distribution 
functions, Vab = V ^aNbVabba- Assume that takes the form of (6), the above equation 
reduces to 

r-Ae 



Wa[l-Za{io)]= Y.VN'b dual,{u)F{u) dio'Re 
^ Jo Jo 

■N{u) + f{-u') , N{u)+f{u') 



X 2 



2 - (w' + Z/)2 



— (oj' — 



(B7) 



and 



Y,^/Nb / dval,{v)F{v) / du;'Re 



V2 - A2(a;') 



X 2 



(a;' + z/)[iV(z/) + /(-a;')] , (a;' - z/)[iVM + /(a;')] 



^2 _ (a;' + z/)2 



dcj'Re 



12 



Aliuj') 



ij^ — ((jj' — Z^)2 



fafe tanh(/3ci;'/2). 



(B8) 



These equations can be solved numerically for the frequency dependence of Z and A. For 
a more transparent demonstration of the physics, we adopt the McMillan approximation [10] 
and look for solutions of the form 



Za{uJ < Aph) = Za{0), Za{uj > Aph) = 1, 

Aaiu < Aph) = A„(0), Aaiu > Aph) = A„(oo). 
This leads to the generalized McMillan formula (7) in main text. 



(B9) 



Appendix C: The phonon signatures 

Conventional signatures of the electron-phonon interaction include the phonon induced 
kink in the normal state dispersion and the phonon shoulder in the tunneling spectra. How- 
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1 2 3 4 5 

(|k|-kF)VF/Aph 



FIG. 9: Solid(dash) line is the renormalized(un-renormalized) electron dispersion, obtained by 
numerical solution of (B7) with parameters Aintra = 0.5 and Aintcr = 0, and the model of FE 
phonon described in Section III. No prominent kink is visible despite significant (factor 1.5) mass 
enhancement at Fermi level. 

ever these features are most pronounced when a'^{u)F{u) have a sharp peak at a character- 
istic phonon frequency. While this is indeed the case for Einstein phonons it is not true for 
the soft phonons under discussion. Here we expect a'^{i')F{u) to have a wide distribution. 
Therefore the above phonon features may not be very obvious. For example using the pa- 
rameters described in the main text a typical renormalized quasiparticle dispersion in the 
normal state near the Fermi surface is shown in Fig. 9. 
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